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ABSTRACT 

We use the integral-field spectrograph SAURON to measure the stellar line-of-sight 
^ ' velocity distribution and absorption line strengths out to four effective radii {Re) in 

OO I the early- type galaxies NGC 3379 and NGC 821. With our newly developed observing 

technique we can now probe these faint regions in galaxies that were previously not 
' accessible with traditional long-slit spectroscopy. We make optimal use of the large 

, field-of-view and high throughput of the spectrograph: by adding the signal of all 

\^ ■ ~1400 lenslets into one spectrum, we obtain sufRcient signal-to-noise in a few hours of 

' observing time to reliably measure the absorption line kinematics and line strengths 

0^ , out to large radius. 

' We find that the line strength gradients previously observed within 1 Rf, remain 

constant out to at least 4 i?e, which puts constraints on the merger histories of these 
galaxies. The stellar halo populations are old and metal-poor. By constructing orbit- 
' based Schwarzschild dynamical models we find that dark matter is necessary to explain 

H I the observed kinematics in NGC 3379 and NGC 821, with 30 - 50 per cent of the total 

matter being dark within 4 R^. The radial anisotropy in our best- fit halo models is 
less than in our models without halo, due to differences in orbital structure. The halo 
also has an effect on the Mg b - Vcsc relation: its slope is steeper when a dark matter 
halo is added to the model. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: individual: NGC 3379, 
NGC 821 — galaxies: kinematics and dynamics — galaxies: haloes — dark matter 



1 INTRODUCTION 

According to the current galaxy formation paradigm, galax- 
ies should be embedded in massive dark matter haloes (e.g. 
Kauffmann & van den Bosch 2002; Springel et al. 2005). For 
spiral galaxies, these haloes have indeed been observation- 
ally confirmed, by determining their influence on the kine- 
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matics of the H I discs that surround most of these galaxies 
(see e.g. van Albada et al. 1985). Early-type galaxies of- 
ten lack these large gas discs, and therefore their dark mat- 
ter haloes are more difficult to probe. Although some early- 
type galaxies do have regular H I structures that have been 
used to identify a dark matter halo (e.g. Bertola et al. 1993; 
Franx, van Gorkom & de Zeeuw 1994; Oosterloo et al. 2002; 
Weijmans et al. 2008) and recent sensitive radio observations 
have shown that H I in these galaxies is more common than 
previously thought (e.g. Morganti et al. 2006), for the ma- 
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jority of early-type galaxies we need to employ other tracers 
to determine the dark matter content. 

Stellar kinematics are a reliable tracer of the gravita- 
tional potential in galaxies. Unfortunately, traditional long- 
slit spectroscopy seldom reaches beyond two half-light or 
effective radii (Re), because the necessary signal-to-noise ra- 
tio (S/N) is difficult to achieve in this regime. In the central 
parts of galaxies this is not an issue, and it has been shown 
by various studies that in early-type galaxies within one Re 
the luminous matter still dominates over the dark matter: 
only 10 - 50 per cent of the total mass is dark (Gerhard et al. 
2001; Cappellari et al. 2006; Thomas et al. 2007). These data 
however do not cover the outer regions of galaxies, where the 
dark matter is supposed to dominate. Models of early-type 
galaxies based on long-slit stellar kinematics extending to 
larger radii are not always conclusive, but indicate that a 
dark halo is required in some individual cases (e.g. CaroUo 
et al. 1995; Rix et al. 1997). 

To probe the gravitational potential at larger radii (out- 
side one -Re), we can resort to gravitational lensing (e.g. 
Koopmans et al. 2006; Gavazzi et al. 2007), X-ray emission 
(e.g. O'SuUivan & Ponman 2004; Humphrey et al. 2006) or 
discrete tracers such as globular clusters (e.g. Cote et al. 
2003; Bridges et al. 2006) or planetary nebulae (e.g. Ro- 
manowsky et al. 2003; Douglas et al. 2007; Napolitano et 
al. 2009). Interestingly, Romanowsky et al. (2003) found 
that models of their sample of three intermediate luminos- 
ity early-type galaxies only required a small amount of dark 
matter, if any at all. Dekel et al. (2005) subsequently ar- 
gued that the observed kinematics of the planetary nebulae 
(PNe) could also be explained by invoking radial anisotropy, 
and therefore do not exclude a dark halo (see also Binney 
& Mamon 1982). Douglas et al. (2007) presented newly re- 
duced PNe kinematics for one of the galaxies of the Ro- 
manowksy et al. sample (NGC 3379) and found the dark 
matter fraction within 5 -Re to be less than 40 per cent, 
which is significantly lower than what is predicted by ACDM 
models. Napolitano et al. (2009) reached similar conclusions 
modeling PNe kinematics in the elliptical galaxy NGC 4494 
and suggested that intermediate-luminosity galaxies are sur- 
rounded by low concentration haloes, while some high-mass 
ellipticals could have more strongly concentrated haloes. De 
Lorenzi et al. (2009) presented mass models of NGC 3379 
based on the PNe kinematics of Douglas et al. (2007), com- 
bined with long-slit and central SAURON integral-field data, 
and concluded that they can fit the data with a range of mass 
distributions, depending on the anisotropy. They did not ex- 
clude a dark matter halo with a strong radially anisotropic 
outer envelope, that is still consistent with merger models 
within the ACDM paradigm. 

We observed two of the galaxies of the Romanowsky 
et al. (2003) sample, NGC 3379 and NGC 821, with the 
SAURON integral-field unit (IFU), to measure the line-of- 
sight velocity distribution (LOSVD) at large radii (3-4 -Re). 
We took advantage of the high throughput and large field- 
of-view of SAURON (41 x 33 arcsec^) and used the IFU as 
a "photon collector": for each observation we co-added all 
spectra within one field to acquire enough S/N to measure 
the LOSVD up to the fourth Gauss-Hermite moment /14. 
These higher order moments are key to breaking the mass- 
anisotropy degeneracy when constructing mass models and 
determining the dark matter content of galaxies (e.g. Ger- 



Parameter 


NGC 821 


NGC 3379 


Morphological Type 


E6? 


El 


Mb (mag) 


-20.74 


-20.57 


Effective B - V (mag) 


0.87 


0.93 


Photometric PA (°) 


32 


68 


Distance modulus (mag) 


31.85 


30.06 


Distance (Mpc) 


23.4 


10.3 


Distance scale (pc arcsec"^) 


114 


50 


Effective radius (arcsec) 


39 


42 


Systemic velocity (km/s) 


1726 


930 



Table 1. Properties of NGC 821 and NGC 3379. Magnitude and 
colour are taken from the Lyon/Meudon Extra galactic Database 
(LEDA). The distance modulus was derived from surface bright- 
ness fluctuation measurements by Tonry et al. (2001). Note that 
0.06 mag is subtracted to adjust to the Cepheid zeropoint of 
Freedman et al. (2001); see Mei et al. (2005), section 3.3, for 
a discussion. The effective radius and position angle are taken 
from Cappellari et al. (2007). 



hard 1993). Apart from the LOSVD, we also obtained line 
strengths from our spectra, allowing us to study line strength 
gradients out to large radii and constraining the proper- 
ties of the stellar halo populations. Both NGC 3379 and 
NGC 821 are well studied galaxies, and some of their prop- 
erties are summarized in Table [1] Their central regions have 
been observed as part of the SAURON survey (de Zeeuw et 
al. 2002). They are classified as fast rotators (Emsellem et 
al. 2007) and display regular stellar kinematics. 

In section 2, we describe our observations and data re- 
duction. We present our spectra and discuss the resulting 
LOSVD in section 3. Section 4 is devoted to measuring line 
strengths and we investigate gradients and the properties of 
the stellar populations at large radii. In section 5 we present 
our mass models of NGC 3379 and NGC 821 (including 
a dark matter halo) and connect the dynamical properties 
with stellar populations by constructing the Mgb - Vcsc re- 
lation. Section 6 summarizes and discusses our results. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed four fields in NGC 3379 with SAURON at the 
William Herschel Telescope at La Palma, Spain, in March 
2005. The fields were centred at 2.6 and 3.5 -Re, on opposite 
sides of the nucleus to allow for an evaluation of the sys- 
temic velocity of the system. They were not placed on the 
major axis of the galaxy but instead on an axis that deviates 
by 23°. This was done to avoid contamination in our outer 
field by the neighbouring galaxy NGC 3384. The positions 
of the fields were chosen such that most of them did not 
contain foreground stars, and they were oriented such that 
the skylenslets of SAURON pointed away from the galaxy 
nucleus, obtaining skyspectra 2 arcminutes away from the 
SAURON field. 

We used a similar strategy to observe six fields in NGC 
821 in September 2006. One field was situated on the minor 
axis, one field at 1 -Re, two fields at opposite sides of the 
galaxy nucleus at 2 -Re and the last two fields on opposite 
sides of the nucleus at 3 -Re . Figure [T] shows the observed 
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Figure 1. Positions of our observed fields in NGC 3379 (left) and NGC 821 (right). The red boxes denote each SAURON field of view. 
The skylenslets (red short thick lines) are aligned with the long side of the SAURON field, at a distance of two arcminutes. SAURON 
was orientated such that the skylenslets pointed away from the galaxy nucleus. The fields are not exactly aligned with the major axis of 
the galaxy (denoted by the dashed line) to minimize the contamination of foreground stars and neighbouring galaxies. The underlying 
images were obtained with the 1.3-m McGraw-Hill Telescope at MDM Observatory. 



Galaxy 


Field 


R/Re 


Tcxp 




S/N 


galaxy contr. 


Comments 










(hr) 


(mag/arcsec^) 




(per cent) 






NGC 3379 


3N 


2.6 


2.0 


23.0 






cloudy conditions, no sif 


;nal 




3S 


2.6 


2.5 


23.0 


29 


15 








4N 


3.5 


7.5 


23.7 


16 


5 


cloudy conditions 






4S 


3.5 


2.0 


23.7 


21 


13 






NGC 821 


minor 


1.0 


1.5 


23.0 


58 


28 








IN 


0.9 


1.5 


21.7 


57 


24 








2N 


1.9 


1.5 


23.3 


12 


7 


cloudy conditions 






28 


1.9 


0.5 


23.3 






cloudy conditions, no sij 


;nal 




3N 


2.9 


4.0 


24.5 






cloudy conditions, no sij 


;nal 




3S 


2.9 


4.0 


24.5 


9 


4 







Table 2. Properties of our observed fields in NGC 3379 and NGC 821. The fields are named according to their distance and orientation 
with respect to the galaxy centre, see Figure [T] The surface brightness fiy is calculated from the MDM images, and can be compared 
to the sky surface brightness, which in V-band is about 21.7 mag/arcsec'^ in optimal conditions at our observed fields (Benn & Ellison 
1998). 



fields in NGC 3379 and NGC 821, and Table [2] summarizes 
their exposure times and exact positions. 

The data were reduced using the dedicated XSAURON 
software (Bacon et al. 2001). The shutter of SAURON mal- 
functioned during the observing run of NGC 821, and did 
not close during read-out of the CCD. This did not affect 
our science frames, whose exposure times of 1800 seconds are 
long compared to the typical read-out time of the CCD (~ 30 
seconds). However, our calibration frames have shorter ex- 
posure times (15 - 60 seconds) and were badly contaminated 
by hght hitting the CCD during read-out. We could there- 
fore not rely on these calibration frames. Instead, we took 
the calibration frames obtained during the observing run of 
NGC 3379 and rotated and translated them so that they 
coincided with the contaminated frames that they would 
replace. The spectra of NGC 821 were then extracted and 
flat-fielded using these frames. 

We checked that the extraction of the spectra was done 
accurately, but we found residuals in the frames after flat- 



fielding. This was to be expected: although the extraction 
mask is rather stable between observing runs, the flat-field 
is not. We therefore constructed a superflat by combining 
our six empty frames with 30 frames of a Lyman a emit- 
ting halo, taken during the same run (see Weijmans et al. 
in prep, for details). Since these 36 frames contain (mostly) 
empty sky and are offset with respect to each other, their 
median, smoothed in the spectral direction, proved to be a 
suitable superflat (see also Bower et al. 2004). Dividing our 
already flat-flelded spectra by this superflat removed most 
of the residuals. Coincidently, we could also construct a su- 
perflat for our NGC 3379 dataset: immediately after our 
observing run SAURON was used in the same setup to ob- 
serve Lyman a emission around high redshift submillimeter 
galaxies. These frames also consist mostly of empty sky, and 
19 of them were used to create a superflat for NGC 3379. 

After superflat-fielding, both datasets were reduced in 
the standard manner, except that we did not subtract the 
sky from the spectra. We will take the contribution of the sky 
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into account when fitting the spectra to extract the stellar 
kinematics (see Section |3}. 



3 STELLAR KINEMATICS 

To extract the stellar kinematics from our spectra, we used 
the penalized pixel fitting method 

(pPXF j3 of Cappellari & 
Emsellem (2004). This method fits a stellar template spec- 
trum, convolved with a line-of-sight velocity distribution 
(LOSVD), to the observed galaxy spectrum in pixel space 
(i.e. logaritmically binned in wavelength). If apart from the 
mean velocity V and velocity dispersion a also the Gauss- 
Hermite moments fts and are fitted, then the solution is 
biased towards a Gaussian. The amount of biasing can be 
controlled, and following the recipe described in Cappellari 
& Emsellem (2004) we found that a penalization factor of 
A = 0.5 gives the best results for our dataset and signal-to- 
noise. 



3.1 Method 

To find a suitable stellar template to fit our spectra, we co- 
added all the spectra of the central SAURON fields of NGC 
821 and NGC 3379, presented in Emsellem et al. (2004). 
We then determined the optimal template for these high 
signal-to-noise spectra as a linear superposition from the 
full MILES library, containing 945 stars (Sanchez-Blazquez 
et al. 2006). For NGC 3379, a template consisting of 19 stars 
was created, and for NGC 821 a template of 21 stars. Tem- 
plate mismatch is negligible, as can be seen in Figure [2] We 
used these templates to measure the stellar kinematics in 
our spectra at large radu in NGC 3379 and NGC 821. Addi- 
tive polynomials were included in our fit to compensate for 
changes in line strengths, which are to be expected over the 
radial range that we cover. We found that varying the stellar 
templates gave similar results for the kinematics within the 
errors. 

Since our spectra are very skydominated (see Table [2]), 
we did not subtract the sky from our spectra, but instead 
included the skyspectra obtained by co-adding the spectra 
observed with the ~170 skylenslets of SAURON (pointing 
~ 2 arcminutes away from the field of view) for each ex- 
posure as separate templates into pPXF. While the stellar 
template was broadened with the LOSVD, at the same time 
the sky spectra were scaled to minimize the x^- The mod- 
eled spectrum Gmod, with x in pixelspace, is then described 
by: 

L N 
Gmodix) = [B® T]{X) + blVl(x) + W„Sr,ix), (1) 

1=0 n=l 

with B{x) the broadening function, T{x) the optimal stellar 
template, 'Pi{x) the L + 1 additive polynomials of degree I 
with coefficients bi, and S„{x) the N sky templates, with 
weight Wn > (see also Eq. 3 in Cappellari & Emsellem 
2004). By fitting the contribution of the sky in this way 
instead of subtracting the skyspectra directly, we allow the 
inclusion of subtle sky variations. We found that especially 

^ available from www-astro. physics. ox. ac.uk/~mxc/idl/#ppxf 
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Figure 2. Co-added spectra of NGC 821 and NGC 3379 (black 
line) and overplotted their optimal stellar templates broadened 
with the best-fitting LOSVD (red line). Residuals (green, shown 
with off'set for presentation purposes) are small, and due to tem- 
plate mismatch. 

for the spectra with a low galaxy contribution (see Table [2]| 
this improved the fits significantly. 

We subtracted the resulting composite skytemplate 
^^^j^ w„Sn{x) from the observed spectrum, and ran pPXF 
again. Errors are estimated by a Monte Carlo method, where 
we measured the kinematics from many realisations of the 
input spectrum, on which we added noise. The noise spec- 
trum was obtained from the difference between our observed 
spectrum and the best-fit modeled spectrum (Gmod). We 
found that varying the amount of subtracted sky by up to 10 
per cent does not affect our measured LOSVD. This shows 
that our LOSVD parameters are robust to errors in our sky 
subtraction. We also checked that the measured LOSVD is 
indeed the LOSVD of the galaxy and not an imprint of the 
solar spectrum from e.g. zodiacal light: when measuring the 
LOSVD of the sky templates themselves we measured a ve- 
locity offset with respect to the galaxy comparable to the 
systemic velocity of the system, and a much smaller broad- 
ening. 

We found that the best fits were obtained by including 
the skyspectra from the skylenslets, taken simultaneously 
with our galaxy spectra. Providing the sky spectra of our 
blank skyfields as templates gave worse results, even though 
these fields covered a larger sky area than the skylenslets and 
therefore had higher S/N. This indicates the importance of 
obtaining simultaneous skyspectra over high S/N, to avoid 
mismatch due to the variability of the night sky. 

Finally, we tested the influence of subtracting residual 
galaxy light that could be present in our sky spectra, since 
the skylenslets were pointing at 6 - 7 -Re in the galaxy. The 
galaxy light at these distances of the nucleus is however 
very faint, approximately 3-4 mag/arcsec^ fainter than in 
the regions where we measure our kinematics. We simulated 
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Figure 3. Spectra of our fields in NGC 3379. The left panels show 
the observed spectra (not skysubtracted) in black, with overplot- 
ted in red the best fit from pPXF (combination of stellar template 
broadened with LOSVD and scaled sky templates). The residuals 
(green) are offset for presentation purposes. Right: skysubtracted 
spectra (black) with the same broadened LOSVD (red) as in the 
left panels. Residuals are again shown in green. 

the effect of subtracting such a weak Gaussian absorption 
line from our observed line and found that the maximal error 
we can introduce in this way is 8 km/s in our measured 
velocity dispersions. This is well within our error bars, and 
we conclude that this effect is negligible. 

3.2 NGC 3379 

It turned out that one of our fields in NGC 3379, 3N, did not 
have sufScient signal to measure the LOSVD. The results for 
the other fields can be found in Figure [3] and Table |31 The 
last four columns for Table|3]show the results for the LOSVD 
if we restrict our fit to the first two moments, instead of 
fitting up to /i4. The results for both fits agree within the 
errors. 

Correcting for barycentric motion, we find that the sys- 
temic velocity Vaya of the galaxy measured from the central 
field (presented by Emsellem et al. 2004) is 930 ± 2 km/s. 
This agrees with V^ya as determined from the 3S and 3N field 
pair: although we could not measure the velocity dispersion 
from the spectrum in 3N, we were able to determine the po- 
sitions of the absorption lines in the spectrum and therefore 
the velocity offset with respect to the 3S field, opposite the 
nucleus. By averaging the two velocities from these field we 
find Kiys = 938 ± 16 km/s. However, repeating this exercise 
for the field pair at 3.5 Re, which are 4S and 4N, we find Vsys 
= 892 ± 23 km/s, which seems marginally deviant from our 
other two measurements, although this is just outside the 
1-a interval. We therefore decided to subtract for each field 
pair their own corresponding systemic velocity. 

In Figure |4] we compare our results with the long-slit 
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Figure 4. Kinematic profiles along the major axis of NGC 3379. 
The small red dots denote the central SAURON data from the 
original survey and the big red dots are our data at large radii. For 
comparison, we also plot the long-slit data of Statler & Smecker- 
Hane (1999) with black stars and the planetary nebulae profiles as 
presented by Coccato et al. (2009) with open diamonds. The hori- 
zontal error bars for the SAURON data at large radii indicate the 
width of the SAURON field-of-view (30 arcsec), while the horizon- 
tal error bars of the PNe data indicate the radial range containing 
68 per cent of the PNe used to calculate the kinematics. 
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Figure 5. Same as FigurefS] but now for our fields in NGC 821. 
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(arcsec) 


(km/s) 
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(km/s) 


(km/s) 


(km/s) 


(km/s) 


3S 


-77 


-77 


8 


14 


155 


14 


-0.04 


0.07 


-0.06 


0.06 


5 


11 


147 


12 


4N 


104 


104 


-19 


19 


126 


34 


0.02 


0.07 


0.05 


0.05 


-16 


19 


134 


35 


4S 


-104 


-104 


22 


14 


124 


25 


0.02 


0.07 


0.02 


0.06 


23 


14 


127 
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Table 3. Stellar kinematics at large radii for NGC 3379. Ax and Ay denote the offsets (East and North, respectively) from the nucleus 
of the galaxy to the centre of our observed fields. The quoted errors are l-tr errors, determined from Monte Carlo simulations, ymean 
and (Tnican are derived by fitting only the first two moments of the LOSVD. 



data of Statler & Smecker-Hane (1999). The radial extent 
of the absorption Une kinematics has been increased by a fac- 
tor two compared to the long-slit data, and our data show 
a smooth continuation of kinematics out to larger radii. We 
see little rotation at large radii, and the dispersion profile 
remains flat out to 3.5 Re, although there is a hint that the 
profile is slightly dropping. We also plot the kinematic pro- 
files of the planetary nebulae, as presented in Coccato et al. 
(2009). We folded these profiles to negative radii, to allow 
comparison with our data at large radii also in these ranges. 
We find a good agreement between the several datasets, al- 
though the rotational velocities found by the planetary neb- 
ulae seem to be a bit higher (~30 km/s) than our findings. 



3.3 NGC 821 

The spectra and results for our fields in NGC 821 can be 
found in Figure [5] and Table |4] Unfortunately, weather con- 
ditions did not allow us to obtain enough signal to measure 
the LOSVD in all our observed fields, but at each radius one 
suitable spectrum could be obtained. The systemic velocity 
from the central kinematics (Emsellem et al. 2004) is 1726 
± 2 km/s (barycentric corrected). This is in agreement with 
Vaya measured from our two field pairs at 1.9 Re (2S and 
2N), where we find a value of 1704 ± 25 km/s. Unfortu- 
nately, we could not measure Vsys from our field pair at 2.9 
Re due to lack of signal in field 3N. However, our field at the 
minor axis of NGC 821 gives an independent determination 
of Vsys, assuming axisymmetry so that here the rotational 
velocity should be zero. We then find Vsys = 1710 ± 7 km/s 
and Vsys ~ 1716 ± 7 km/s, depending on whether we fit up 
the LOSVD up to or only the first two moments, respec- 
tively. Since all these measurements of Vsys agree within the 
errors, we fixed Vsys such that the velocity measured on the 
minor axis is zero. 

Figure [6] compares our result with the long-slit data 
of Forestell & Gebhardt (2008). We note that especially in 
the centre there seem to be some deviations between this 
dataset and the SAURON data. Differences in sampling be- 
tween the two datasets (nearly-round bins for the SAURON 
data versus long, elongated bins for the long slit data) could 
in principle result in differently-weighted values, but we ex- 
pect that this effect is small. The most likely cause for the 
observed deviations are differences in fitted wavelength in- 
terval and stellar templates. Also shown are the results from 
Coccato et al. (2009), who re-analysed the planetary nebu- 
lae data that were previously presented by Romanowsky et 
al. (2003). These data are in good agreement with our stellar 
kinematics at large radii. 



4 LINE STRENGTHS 

To measure the absorption line strengths in our spectra, we 
follow the procedure outlined in Kuntschner et al. (2006). 
Briefly, we can measure three Lick indices in the SAURON 
wavelength range: H/9, Fe5015 and Mg6, deflned in Table 2 
of Kuntschner et al. (2006). We do not need to remove emis- 
sion lines from our spectra, since our galaxies show no sign 
of ionised gas outside the central regions (Sarzi et al. 2006). 
To calibrate to the Lick/IDS system, we broaden our spectra 
to an instrumental resolution of a = 212 km s~^. We cor- 
rect our measured line strengths for the broadening caused 
by the LOSVD and we apply the offsets from Table 3 in 
Kuntschner et al. (2006) to correct for differences in the 
shape of the continuum, since the Lick/IDS spectra have 
not been flux calibrated. 

We noticed that the continuum shape blue-ward of the 
H/9 feature in our spectra was affected by bad flatfielding or 
wrong sky subtraction. We therefore applied the same con- 
tinuum correction as was done in Kuntschner et al. (2006) 
and fitted an 11**^ order multiplicative polynomial together 
with an optimal stellar template to our spectra. The tem- 
plate had been broadened with the observed LOSVD. We 
then divided our observed spectrum by this polynomial, be- 
fore determining the H/3 index. 

4.1 Results 

We present our resuffs for NGC 3379 and NGC 821 in Ta- 
ble [S] Errors are estimated with Monte Carlo simulations, 
where we add noise to our spectra (defined by the difference 
between our observed and fitted spectra, see Figures |3] and 
[S]) and vary the bandpasses of the index and the continuum 
to take the uncertainty in radial velocity into account. 

A possible source of contamination for our line strength 
measurements is the [O ill] emission of planetary nebulae. 
From the catalog of Douglas et al. (2007), each of our fields in 
NGC 3379 contains one to two planetary nebulae. Following 
Sarzi et al. (2005), their Equation 1, we should be sensitive 
to emission lines with EW > 0.7 A. We carefully checked 
the individual spectra in each cube, prior to co-adding, and 
found no evidence for such emission. 

We also investigated the effect of sky subtraction on the 
determinations of line strength. By varying the amount of 
sky that we subtracted from our spectra, we found that our 
measurements of H/3 and Fe5015 are quite robust against an 
error in sky subtraction, even if the skylevels are varied by 10 
per cent. Mg6 is however already severely affected by a sky 
subtraction error of only two per cent (see Figure (Tj . This 
may be due to the proximity of the solar Mg b absorption 
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Table 4. Stellar kinematics at large radii for NGC 821, similar to Table [3] 
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Figure 6. Kinematic profiles along the major (left) and minor (right) axis of NGC 821. As in Figure |4] small red dots are the central 
SAURON data of the original survey, big red dots are our data at large radii and open diamonds denote planetary nebulae kinematics 
from Coccato et al. (2009). The black stars are long-slit stellar kinematics from Forestell & Gebhardt (2008). Horizontal error bars are 
as described in Figure |4] 
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Table 5. Line strengths indices in our observed fields in NGC 
821 and NGC 3379. 



feature present in the sky spectrum. Our Fe5015 and Mgfe 
indices at 2 and 3 Re in NGC 821 are unrealistically high, 
and therefore could be suffering from this problem, although 
variations in the continuum shape of the spectrum also play 
a role (see next section). The indices at 1 in this galaxy 



are less affected, and are in agreement with Proctor et al. 
(2005). They determined line strengths from long-slit data 
in this galaxy out to 1 i?e . We are not aware of studies in the 
literature where line strengths have been determined outside 
1 Re for these galaxies, to compare our results with. 

4.2 Line strength gradients 

Line strength gradients are key instruments to study the 
formation history of galaxies, since different formation mod- 
els predict different gradients. Monolithic collapse models 
(Carlberg 1984) produce steep metallicity gradients and a 
metal-rich nucleus, while hierarchical models show shallower 
gradients, as mergers dilute existing line strength gradients 
(White 1980). However, subsequent starformation induced 
by gas-rich mergers can steepen the gradients again, at least 
in the central regions (Hopkins et al. 2009a). 

Our measurements allow us to study line strength gra- 
dients out to large radii. Many previous studies of galaxy 
samples (e.g. Davies, Sadler & Peletier 1993; Kuntschner et 
al. 2006; Sanchez-Blazquez et al. 2007) have been limited 
to ~ 1 Re, but showed that typically for galaxies without 
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Figure 7. The effect of sky subtraction on our line strength measurements, for the field 3S in NGC 3379 (left) and the IN field in 
NGC 821 (right). The x-axis denotes the fraction of sky subtracted. The dashed line indicates the line index when the sky is totally 
removed, while the grey box denotes the l-a interval. H/3 and Fe5015 are robust against an error in sky subtraction up to 5 per cent, 
though oversubtracting is more severe than undersubtracting. Mg b is very sensitive to sky subtraction errors. 
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Figure 8. Line strength indices H/3, Fe5015 and Mgb (in magnitudes) as a function of radius. Left panel: NGC 3379, right panel: NGC 
821. Black triangles show indices measured on the SAURON central field, taken from a slit positioned on the major axis of the galaxy, 
with a width of 2 arcsec. Red filled dots show our large radii results, and the horizontal error bars denote the size of the SAURON field of 
view. The points at 2 and 3 Re in NGC 821 cannot be considered reliable, probably because of continuum shape problems, and therefore 
are not shown here. The blue open squares denote our measurements directly from templates of single stellar population models, see text 
for details. Overplotted is the best fitting linear relation between index and log(radius), determined within Re /A. 



contributions from young stellar populations, H/3 remains 
flat or increases slightly with radius, while the metal line 
strengths have negative gradients. 

To calculate the line strength gradients in our two galax- 
ies, we follow Kuntschner et al. (2006) and first convert our 
indices from A to magnitudes: 



index' = -2.5 log(l-^), (2) 



with AA the width of the index bandpass in A. Indices ex- 
pressed in magnitudes will from now on be indicated by a 
prime sign [']. Line strength gradients are then defined as: 
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Aindex' — 



Sindex' 
5log{R/R, 



(3) 



We show in Figure |S] our line strength indices in magnitudes 
as a function of radius, with red dots. Also plotted are the 
values for the central SAURON fields (Kuntschner et al. 
2006), which fall within a 2 arcsec wide slit positioned along 
the major axis of the galaxy. We fitted a linear relation be- 
tween each index and log(radius) within i?e/4, to determine 
the gradient. These gradients are overplotted in the figure, 
and are compared to our measurements at large radii. 

With our new data, we confirm in NGC 3379 the trends 
mentioned above, that H/3 mildly increases and Fe5015 and 
Mg6 decrease with radius. This indicates that the metallic- 
ity in the galaxy keeps decreasing linearly with log(radius) 
when going from the centre out to larger radii. The gra- 
dients for H/3 and Mgfo seem to remain constant out to 4 
Re, but Fe5015 falls off more steeply. This would imply an 
unrealistically high a-enhancement at large radii. To check 
whether our Fe5015 measurements are affected by contin- 
uum shape problems, we did the following test. We remea- 
sured the kinematics in our spectra with pPXF, but this 
time used the stellar model library of Vazdekis et al. (2007) 
based on the MILES library (-1.68 < Z/H < 0.0; 1 Gyr 
< age < 17.78 Gyr) and multiplicative instead of additive 
polynomials, to allow for a varying continuum shape. The 
resulting kinematics were consistent with our previous find- 
ing. The line strengths were then measured on the compos- 
ite best-fitting stellar template, without the multiplicative 
polynomials. The results are also shown in Figure |8] as blue 
squares. For H/3 and Mg6 both the direct measurements and 
the measurements on the stellar templates give similar re- 
sults, but the values for Fe5015 are higher and consistent 
with the observed gradient in the central part of the galaxy. 
We therefore conclude that indeed our measurements for 
Fe5015 suffer from changes in continuum shape, but that 
our values for H/3 and Mg6 are robust. For NGC 821 we see 
a similar effect. 

Our observation that the line strength gradients are 
constant out to large radii provides constraints for the 
merger and star formation history of the galaxy. After a 
gas-rich merger, star formation in the central regions of the 
remnant is expected to steepen the gradients, while at larger 
radii violent relaxation fiattens them, though over time (~ 3 
Gyr) the steep gradients in the central part get weaker (Hop- 
kins et al. 2009a; Hopkins et al. 2009b). 



model predictions from Thomas et al. (2003) is zero. In Fig- 
ure |9] we plot this index versus H/3 and compare with the 
models. Even though the models of Thomas et al. (2003) and 
Schiavon (2007) are constructed independently, using differ- 
ent stellar libraries and fitting functions, they give the same 
results for our data at large radii. For NGC 3379 we find 
that at 3 - 4 _Re, the stellar population is consistent with an 
old (12 Gyr) population, and metal-poor, with [Z/H] slightly 
below 20 per cent of the solar metallicity. We note however 
that the uncertainty in these values is rather large. For NGC 
821, we find that at 1 Re the stellar population is of the same 
age and metallicity range as obtained from the outer bins of 
the SAURON central field (~ 0.6 Re). 

Stellar population models predict a decrease in stellar 
mass-to- light ratio M,/L if the metallicity of the stellar pop- 
ulation decreases, since the stars then become bluer and 
therefore brighter in the optical. For instance, for a change 
in metallicity from 0.0 to -0.84 at a constant age of 10 Gyr 
(consistent with our observervations), M*/L decreases by 
about 23 per cent in the models of Maraston (2005). How- 
ever, for increasing stellar age at constant metallicity, M, /L 
also increases. Given this degeneracy and uncertainties, we 
adopt a constant stellar M/L while constructing the dy- 
namical mass models that we present in the next section. 
Furthermore, M^/L depends strongly on H/3 (see e.g. Cap- 
pellari et al. 2006). Since the H/3 profiles in our galaxies are 
nearly fiat out to large radii, variations in M«/L are most 
likely small. 



5 DYNAMICAL MODELS 

To explore whether our data of NGC 3379 and NGC 821 are 
consistent with a dark matter halo, we model these galaxies 
with the triaxial Schwarzschild code presented by van den 
Bosch et al. (2008). Orbits are calculated within an a pri- 
ori specified triaxial potential, and a superposition of orbits 
is determined that best fits the observed kinematics, using 
the observed surface brightness and self-consistency as con- 
straints. This way, the viewing directions, stellar mass-to- 
light ratio M,/L, central black hole mass M, as well as the 
internal orbital structure can be determined (see also van de 
Ven, de Zeeuw & van den Bosch 2008; van den Bosch & van 
de Ven 2008). Within the limits of the observed photometry 
and LOSVD, the velocity anisotropy of the galaxy can vary 
freely. 



4.3 Stellar population models 

Using the stellar population models of Thomas, Maraston 
& Bender (2003) and Schiavon (2007), we explore the sin- 
gle stellar population (SSP) equivalent age and metallicity 
in NGC 3379 and NGC 821 at large radii. Following the 
aproach of Thomas et al. (2003) we define an abundance 
ratio insensitive metallicity index similar to their [MgFe]' 
index by using the Mg6 and Fe5015 indices available from 
SAURON. [MgFeSO]' is defined as: 



[MgFeSO]' = 0-69xMg6 + Fe5015^ 



(4) 



The scaling factor for the Mgfe index was optimized such 
that the mean difference between solar and non-solar ratio 



5.1 Including a dark halo 

We extended the triaxial Schwarzschild code by adding a 
dark matter halo to the gravitational potential, which al- 
ready consisted of a stellar contribution inferred from the 
observed photometry and a Plummer potential for the con- 
tribution of the central black hole (see van den Bosch et al. 
2008 for more details). We parametrise the dark halo with 
a spherical NFW profile (Navarro, Frenk & White 1996) ; 



p{r) 



(5) 



with ps the characteristic density of the halo and a char- 
acteristic radius. The potential generated by this density 
distribution is then given by: 
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Figure 9. H/3 index against [MgFe50]' (both in A) in NGC 3379 (top) and NGC 821 (bottom). Loft panels show the Thomas et al. (2003) 
stellar population models (solid and dashed lines) and right panels the Schiavon (2007) models. Black dots indicate measurements from 
the SAURON central field, while the coloured dots are averaged along isophotes (see inset for colour coding). The black filled squares 
show the data at large radii, revealing an old (12 Gyr) and metal-poor (below 20 per cent solar metallicity) population for both models, 
though the uncertainties in these values are large. 



$(r) = -47rGpsr,'— ln(l + — ), (6) 
r rs 

where G is the gravitational constant. 

We are interested in the presence and acceptable mass 
range of the halo; modeling the shape of the halo in detail 
is beyond the scope of this paper. We therefore fix the con- 
centration c = r2oo/rs of the halo to 10, as indicated by 
cosmological simulations (Bullock et al. 2001). Here r2oo is 
the radius of the halo within which the mean density has 
dropped to 200 times the critical density pcrit- It can be 
shown that 

Ps _ 200 

Pcrit 3 ln(l + c)- c/(l + c) ' ^ ' 

such that c determines ps in the potential of Equation [S] We 
vary the halo mass, 

M200 = 47rp,r? f ln(l + c) - -^1 , (8) 
'- 1 + C-' 



which is the enclosed mass within r2oo, yielding Vs- We only 
fit for M200 and keep M^/L fixed to the maximum value al- 
lowed by the data, therefore fitting a minimal halo or max- 
imal spheroicJJ. The influence of the black hole is negligible 
at large radii and therefore we also do not fit for M, but 
keep its value fixed. 

5.2 Modelling NGC 3379 and NGC 821 

NGC 3379 has a small kinematic misalignment of 5 ± 3° 
(e.g. Statler & Smecker-Hane 1999) and shows signs of a 
kinematic twist (Krajnovic et al. 2008), which indicate that 
the galaxy is not perfectly axisymmetric (see also Capaccioli 
et al. 1991; Statler 2001). Van den Bosch & de Zeeuw (2009) 

^ Comparable to the maximal disc model that is often invoked 
when modeling dark haloes of spiral galaxies 
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therefore constructed triaxial Schwarzschild models of this 
galaxy. Their best-fitting triaxial model (axis ratios p = 0.95, 
q = 0.81) is nearly face-on, with an inclination of ~ 40°. 
This model is nearly spherical in the centre and becomes 
close to oblate at large radii, which is in agreement with 
Statler (2001) and De Lorenzi et al. (2009). 

We adopt the best-fit model of van den Bosch & de 
Zeeuw (2009) for our halo modeling. The model has a black 
hole mass of M, = 4.0 x 10* M© and a dynamical mass- 
to- light ratio M/Li — 2.9 Mq/Lj^q, which we adopt as 
the maximally allowed stellar mass-to-light ratio M,/L. The 
stellar potential is represented by a Multi- Gaussian Expan- 
sion (MGE, Emsellem, Monnet & Bacon 1994; Cappellari 
2002) model based upon /-band imaging], presented by 
Cappellari et al. (2006). We use the same central SAURON 
kinematics (up to /ig) extracted with the MILES library (see 
Shapiro et al. 2006; Cappellari et al. 2007; van den Bosch & 
de Zeeuw 2009) to cover the inner Re of NCC 3379 and add 
our measure points at large radii. For extra spatial coverage, 
we also include the symmetrized data of the four long-slit ob- 
servations of Statler & Smecker-Hane (1999) (see Figure fTO)) . 
The large mosaic of central SAURON data provides many 
more kinematic constraints in the central region than the 
four slits. We therefore exclude the inner 20 arcsec of each 
slit from the fit. The planetary nebulae measurements are 
not based on integrated stellar light, as our SAURON and 
long-slit data, and therefore for consistency not included. 

We use a similar procedure to construct mass models 
of NCC 821. This galaxy is a very flattened system (E6) 
and shows no signs of non-axisymmetry in its central part, 
even though planetary nebulae kinematics suggest twists at 
larger radii (Coccato et al. 2009). Our dataset however pro- 
vides no constraints on the behaviour of the velocity fleld at 
these scales. We therefore deproject this galaxy axisymmet- 
rically, assuming an edge-on view and use the triaxial code in 
the axisymmetric limit. We adopt the edge-on MGE model 
of NCC 821 based upon ground-based ii-band photome- 
trjQand space-based imaging with HST/WFPC2 in F814W- 
band, as in McDermid (2002) and Cappellari et al. (2006). 
We first construct a model without a dark halo, based only 
on kinematics (up to fee) of the central SAURON dataset, 
extracted with the MILES library (Cappellari et al. 2007). 
Our best fit has a stellar mass-to-light ratio of 2.5 Mq/Lj^q 
with a black hole mass of M, = 2.1 x 10** from the M.- 
cr relation. We then add a dark NEW halo to this model, 
with a fixed concentration of c = 10. The SAURON cen- 
tral dataset is used to cover the central part of the galaxy, 
while our four measurepoints at large radii are included to 
cover the outerparts. To avoid a possible conflict between the 
central SAURON data and the long-slit data of Forestell & 
Gebhardt (2008), we decided to not include the latter in our 
model, as there seem to be some deviations between the two 
datasets (see Figure [HJ. 

Our best fitting model for NCC 3379 requires a halo 
mass of Af200 ~ 1-0 x W^^ Mq, as compared to a total stellar 
mass of 1.0 x 10^^ M©. In Figure [T^ we show the kinematics 



^ Obtained with the 1.3-m McGraw-Gill Telescope of the MDM 
Observatory, Kitt Peak (Falcon-Barroso et al. in prep) 
^ Obtained with the 0.9-m telescope at the Cerro Tololo Inter- 
American Observatory (CTIO), Chile 




Figure 10. Overview of the datasets included in our dynamical 
model of NGC 3379. The dotted lines denote isophotes of the 
galaxy at 1, 2, 3 and 4 Re. The (red) box in the centre denotes 
the central SAURON field (Emsellem et al. 2004), while the other 
boxes show the positions of our data at large radii. The (blue) 
solid lines show the orientation and extent of the long-slit data 
(Statler & Smecker-Hane 1999). 

of our best-fit model compared to our data at large radii 
and the long-slit kinematics on the major axis. The fit is 
excellent, even for the higher order moments hs and /14. A 
model without dark halo clearly underpredicts the observed 
velocity dispersion and hi outside 1 Re - For comparison also 
a model with a halo heavier than our best fit model is shown: 
here the velocity dispersion is overestimated. In Figure [TT] 
we show the values of our best fitting model. We see that 
with 3-cr confidence, M200 = l-0io;3 ^ 10^^ Mq under the 
assumptions mentioned above, while a model without dark 
halo is strongly excluded. 

We obtain similar results for NGC 821, where we also 
exclude the absence of a dark halo (see Figure [T5|l. The halo 
mass in our best-fit model is M200 ~ 6.514 3 ^ IO'^^Mq, 
while the total stellar mass is 1.1 x 10^ '^Mq. The kinematics 
of the fit compared to the observed kinematics are shown 
in Figure [T4l where for comparison also a model without a 
dark halo and a model with a too massive halo are shown. 

We conclude that dark matter is present both in 
NGC 3379 and NGC 821. Even within 1 Re in NGC 3379, 
8 per cent of the total mass is dark, while at 4 Re this frac- 
tion has increased to 34 per cent. NGC 821 is more dark 
matter dominated, with a dark matter fraction of already 
18 per cent within 1 Re, increasing to 49 per cent within 
4 Re- Note however that these values are based on mini- 
mal halo assumptions, and therefore should be seen as lower 
limits. Also, varying halo shape and profile can cause these 
fractions to change, they therefore should be strictly seen as 
indicative. 

5.3 Orbital structure 

In the top panels of Figure [15] we show the orbital anisotropy 
in our dynamical models of NGC 3379 and NGC 821. We 
characterize the orbital anisotropy as 



with (r, 8, (f)) standard spherical coordinates. Our best-fit 
model of NGC 3379 is nearly isotropic within 2 Re, and 
at larger radii the radial anisotropy increases. Our best-fit 
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Figure 11. Ax^ levels for the Schwarzschild models we ran for 
NGC 3379. The only free parameter in these models was A/200 
(total halo mass). The model without dark halo is relocated at 
3 X W^^Mq, to place it within the plot. 
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Figure 12. Best-fit model for NGC 3379 (bold Une, M200 = 1.0 X 
lO^^Afo) compared to our datapoints at large radii (red filled 
circles) and long-slit data on the major axis (black stars). Also 
indicated are a model without halo (dashed line) and a model with 
a heavy halo (dashed-dotted line, M200 = 5.2 X IO^^Mq). These 
models respectively underpredict and overpredict the observed 
velocity dispersion. 



model of NGC 821 is already mildly radially anisotropic in 
the central regions, but also shows a small increase with ra- 
dius. Cappellari et al. (2007) presented orbital anisotropy 
values for both NGC 3379 and NGC 821 from axisymmetric 
models, based on the same central SAURON data and MGE 
models as we use in this paper. Our central anisotropy val- 
ues are consistent with theirs for both galaxies. Since Cap- 
pellari et al. (2007) used a different, independent modeling 
code, this agreement gives confidence in the robustness of 
the anisotropy parameters presented here. 
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Figure 13. Same as Figure [TTI but now for NGC 821. 
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Figure 14. Best-fit model for NGC 821 (bold line, Af2oo = 
6.5 X 10^^ Mq) compared to our datapoints at large radii (red 
filled circles) and the central SAURON data (open triangles). The 
dashed line indicates a model without halo and the dashed-dotted 
line a model with a too heavy halo (A^200 = 5.5 X IO^^Mq). 



Within 1 Re, the orbital anisotropy changes only 
marginally between our models with and without halo, con- 
sistent with the findings of Thomas et al. (2009). At larger 
radii the radial anisotropy increases, and the increase is 
smaller for models with heavier haloes. This is contrary to 
what is seen in spherical models. In order to fit the observed 
mildly declining or constant velocity dispersion profile with- 
out invoking extra dark matter at large radii, the orbits have 
to be tangentially anisotropic, while adding a dark halo re- 
sults in more radially anisotropic orbits (e.g. Richstone & 
Tremaine 1984; Dekel et al. 2005; de Lorenzi et al. 2009). 

The detailed kinematics in the inner part of the galaxy 
together with our observed positive in the outer parts 
already tightly constrain the fraction of box and tube or- 
bits in our triaxial potential, and therefore the anisotropy. 
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Figure 15. Orbital structure in NGC 3379 (left) and NGC 821 (right). Top panel: velocity anisotropy /3r as a function of radius (see 
Equation (Ojl . The lines are as in Figures 1121 and 1141 with the solid line for our best- fit model, the dashed line for the model without halo 
and the dashed-dotted line for the model with the too heavy halo. Bottom panel: fraction of orbit types as a function of radius for the 
models of the top panel. Blue lines denote the fractions of short axis tubes, red the fractions of long axis tubes and black the fractions 
of box orbits. 



Box orbits have no net rotation, and contribute to the ra- 
dial anisotropy of the system, requiring /i4 > 0. For tube 
orbits to reproduce the observed nearly zero mean rotation 
at large radii, they have to be counter-rotating, increasing 
the velocity dispersion. However, they would at the same 
time increase the tangential anisotropy and therefore yield 
hi < 0. Since we observe a positive /i4 at large radii, box 
orbits are expected to contribute significantly to the orbital 
structure in our models without dark matter, as confirmed 
in Figure [151 As a result, the models without dark halo are 
strongly radial anisotropic at large radii. Adding a spheri- 
cal dark matter halo to the system makes the gravitational 
potential less triaxial. This will decrease the phase space 
volume of box orbits and naturally results in a larger popu- 
lation of tube orbits (Hunter & de Zeeuw 1992). The radial 
anisotropy of the system will decrease, as is illustrated in 
Figure 1151 This decrease is more pronounced in NGC 821 
than in NGC 3379, since our best-fit model of NGC 821 is 
more dark matter dominated than our NGC 3379 model. 

5.4 Mgfe - Vcsc relation 

We here investigate the connection between the stellar pop- 
ulations and the mass distribution of the galaxy. Franx & 
Illingworth (1990) showed that there is a tight relation be- 
tween the local colour (tracing metallicity) and the local 
escape velocity Vcsc (tracing the gravitational potential $) 
in early-type galaxies, with Vcsc given by: 

Vcsc = V2i*l- (10) 

This relation was confirmed by Davies et al. (1993), who 
used Mg 2 as a tracer for metallicity. Scott et al. (2009) find 
correlations between line strength indices and Vesc for the 
early-type galaxies in the SAURON sample, based on mea- 
surements within ~1 f?e. They also show that for Mg6 this 
relation does not only hold locally in one galaxy, but is also 
a global relation between different galaxies. 

With our data at large radii we can now investigate 
the behaviour of this relation out to 4 Re- We determine 
for our best-fit model of NGC 3379 the contribution to the 



gravitational potential of the dark halo using Equation |6l 
The contribution of the stellar mass is inferred from the 
MGE model, following the formularium in the appendix of 
Cappellari et al. (2002) and the black hole is included as 
a Plummer potential. We plot our result in Figure [T51 and 
find that the observed relation within 1 Re continues out to 
at least 4 Rc . The influence of the dark halo on the Mg b - 
Vcsc relation becomes clear when comparing our results to a 
model without a dark halo, where the only contribution to 
the potential comes from the MGE model and the black hole 
(open symbols in Figure I16p . Not only is there an offset in 
Vcsc with respect to the model with the dark halo, but also 
the slope changes, an effect which was also noted by Franx & 
Illingworth (1990). The addition of extra mass at larger radii 
by the dark halo deepens the potential well, which means 
that even particles in the centre of the galaxy will need a 
larger Vase to climb out of this well, causing a shift in Vcsc 
at all radii. However, particles in the outskirts of the galaxy 
are relatively more affected by the addition of the dark halo 
than particles in the more central regions: if the halo had not 
been present, escaping from the potential well at large radii 
would have been relatively easy. Therefore, the shift in Vcsc 
is larger at larger radii, something that is also illustrated by 
Figure 1161 The net result is a change in slope. We find that 
the slope of the Mg 6 - Vcsc relation in NGC 3379 changes in 
logscale from 0.31 to 0.46, when taking the dark halo into 
account. If instead of considering the infinite halo mass of 
the NEW profile we only take the contribution of the halo 
within r2oo into account when calculating the gravitational 
potential, the slope becomes 0.43 in logscale. 

We note that in both the models with and without halo, 
there is a local relation between Mg6 and Vcsc- However, 
when looking at a larger sample of galaxies and studying 
the global Mgb - Vcsc relation the contribution of the halo to 
the potential should be taken into account to interpret the 
slopes, since the relative contribution of the dark matter 
to the potential compared to the luminous one may vary 
from galaxy to galaxy, e.g. depending on galaxy mass or 
environment. 
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Figure 16. Mgb (in A) versus local escape velocity Vesc (km/s) in 
NGC 3379. Black filled triangles denote points from the central 
SAURON data, while the black filled dots are our large radii 
data. Overplotted is the gradient fitted to the central dataset. 
With open symbols (in blue) we show the Mg6 - Vcsc relation 
calculated from our model without a halo. This relation is offset 
from our best-fit model and also has a different slope, illustrating 
the influence of the dark matter. 



6 SUMMARY AND CONCLUSION 

We used the SAURON IFU as a "photon collector" to obtain 
stellar absorption line kinematics and line strengths out to 
four efltective radii in the early-type galaxies NGC 3379 and 
NGC 821. By co-adding the individual spectra of each lenslet 
into one spectrum for each observed field, we were able to ob- 
tain sufficient signal-to-noise in the outskirts of these galax- 
ies, which are too faint to observe with traditional long-slit 
spectroscopy. The stellar velocity profiles are measured up 
to the fourth Gauss-Hermite moment {hi), which allows us 
to break the mass-anisotropy degeneracy when constructing 
dynamical models of our galaxies. The kinematics show a 
smooth continuation of the profiles measured within 1-2 Re 
with SAURON and long-slit spectroscopy, and are consis- 
tent with planetary nebulae kinematics. The line strengths 
also show a continuation of the gradients observed within 
1 although our Fe5015 measurements probably suffer 
from variations in the continuum shape. The stellar halo 
population is old and metal-poor. The continuation of the 
line strength gradients and the Mg 6 - Vcsc relation out to at 
least 4 Re provides constraints for the star formation and 
merger history of early-type galaxies (e.g. Hopkins et al. 
2009a; Hopkins et al. 2009b; Spolaor et al. 2009). 

Constructing dynamical orbit-based models, we show 
that even assuming a maximal stellar contribution, both 
NGC 3379 and NGC 821 require a dark halo to fit the 
observed kinematics. We parametrised the halo contribu- 
tion with a spherical NFW profile, fixing the concentration 
to c = 10. We therefore cannot give a definite number for 
the precise mass content, shape and anisotropy of the halo. 
Varying c in our halo profile does not have a large effect on 
the fitted dark matter content, as an increase of c will be ac- 
companied by a decrease in M200, as the dynamical models 
constrain the density inside the region for which kinematics 
are available. Only with data at even larger radii (preferably 
out to ~ 100 Re, accessible with e.g. weak gravitational lens- 
ing or satellites) we would be able to break the degeneracy 



between c and M2oo- In that sense, our results should be 
viewed as an estimate of the minimal contribution of the 
dark halo to the total mass budget in the galaxy, and we 
then see that even at small radii (within 1 Re), the dark 
matter is already responsible for 10 - 20 per cent of the 
total mass. These percentages will increase if we relax the 
minimal halo assumption, but we then need reliable mea- 
surements of the stellar mass-to-light ratio and its poten- 
tial variation with radius. With better spatial coverage of 
the galaxy, either by mosaicing or by employing larger IFUs 
(such as VIRUS-P and in the future MUSE) and combining 
the so-obtained stellar kinematics with other datasets trac- 
ing dark matter at different regimes (e.g. H I discs or rings, 
hot gas, gravitational lensing and/or discrete tracers such as 
planetary nebulae) we will be able to build more sophisti- 
cated dynamical models, to explore the shape and mass of 
the dark haloes around early-type galaxies. 
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